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STABILITY  OF  MOTIONS  OF  THERMOELASTIC  FLUIDS 


C.  M.  Dafermos 

V 

ABSTRACT : It  is  shown  that  the  second  law  of  thermodynamics  in- 
duces uniqueness  and  continuous  dependence  upon  initial  state  and 
supply  terms  of  smooth  thermodynamic  processes  of  thermoelastic 
fluids  within  the  broader  class  of  thermodynamic  processes  with 
shock  waves. 
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1.  Introduction 


The  intent  of  this  article  is  to  contribute  to  the  project  of 
elucidating  the  relationship  between  the  second  law  of  thermodynamics 
and  stability.  The  early  work  of  Duhem  [1]  and  its  subsequent 
developments  by  Ericksen  [2,3],  Coleman  and  Dill  [4]  and  Gurtin  [5] 
indicate  that  the  second  law  (the  Clausius-Duhem  inequality)  induces 
Liapunov  stability  of  the  equilibrium  process  in  thermomechanics. 

In  a different  direction,  the  author  [6]  shows  that  the  second  law 
renders  uniqueness  of  thermoelastic  processes  as  well  as  continuous 
dependence  upon  the  initial  state  and  the  supply  terms. 

Here  we  continue  the  investigation  in  [6]  by  examining  the 
implications  of  the  second  law  on  the  stability  of  adiabatic  processes 
of  thermoelastic  fluids.  This  problem  is  of  some  interest  for  the 
following  reason:  The  nonlinearity  and  hyperbol icity  of  the  field 
equations,  expressing  the  balance  laws  of  mass,  momentum  and  energy 
in  the  theory  of  thermoelastic  fluids,  causes  the  breakdown  of 
smooth  solutions  and  the  development  of  shock  waves.  Thus,  the 
class  of  smooth  functions  is  far  too  narrow  to  encompass  all 
processes  of  physical  interest.  Careful  mathematical  investigation 
[7]  reveals  that  processes  should  be  defined  in  the  set  of  func- 
tions of  bounded  variation,  in  the  sense  of  Tonelli  and  Cesari, 
that  is  the  class  of  velocity,  density  and  entropy  fields  whose 
gradients  are  Borel  measures.  The  reader  who  does  not  care  to  get 
involved  with  technicalities,  may,  without  much  loss,  visualize 
these  processes  as  smooth  except  on  a family  of  propagating  surfaces 
(shocks)  across  which  the  velocity,  density  and  entropy  fields 
experience  jump  discontinuities. 


Within  the  class  of  processes  of  bounded  variation,  the  initial 
value  problem  admits,  in  general,  many  solutions,  i.e.,  there  may 
be  many  processes  that  originate  from  the  same  state  and  satisfy 
the  balance  laws  of  mass,  momentum  and  energy  for  the  same  body 
force  and  heat  source.  Interpreted  as  an  admissibility  criterion, 
the  second  law  of  thermodynamics,  in  the  form  of  the  Clausius-Duhem 
inequality,  rules  out  some  but  not  necessarily  all  extraneous 
processes.  Consequently,  in  order  to  single  out  the  physically 
relevant  process  one  has  to  employ  stronger  admissibility  criteria 
[8],  In  this  paper,  however,  we  show  that,  somewhat  surprisingly, 
whenever  a smooth  process  exists  then  it  is  unique  and  stable  within 
the  broader  class  of  processes  of  bounded  variation  that  satisfy  the 
Clausius-Duhem  inequality.  In  other  words,  for  as  long  as  one  is 
dealing  with  smooth  processes,  the  second  law  of  thermodynamics  in 
its  traditional  form  is  sufficiently  powerful  to  rule  out  all 
extraneous  processes. 

The  proof  of  the  aforementioned  result  is  based  on  an  energy 
inequality,  derived  in  Section  3,  in  the  spirit  of  the  "entropy" 
estimate  established  by  DiPerna  in  his  interesting  paper  [9]  on 
uniqueness  of  solutions  of  quasilinear  hyperbolic  systems. 
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2.  Thermoelastic  Fluids 


We  begin  with  a review  of  the  classical  theory  of  thermo- 
elastic fluids.  We  will  be  concerned  here  exclusively  with 
adiabatic  processes  so  that  the  state  variables  are  velocity  v, 
specific  volume  t (the  inverse  of  density  p,  i.e.,  i = 1/p) , and 
entropy  n.  In  terms  of  these,  internal  energy  e,  pressure  p, 
and  temperature  0 are  determined  via  constitutive  equations: 

e = e*(T,n),  P = p*(t,n),  0 = e*(T,n)  (2.1) 


P*(T,n) 


e*(*,n) 


(2.2) 


An  admissible  thermodynamic  process  is  determined  by  velocity, 
specific  volume  and  entropy  fields  (v(x,t),  x(x,t),  n(x,t))  that 
satisfy  the  balance  laws  of  mass,  momentum  and  energy  as  well  as  the 
second  law  of  thermodynamics,  viz., 
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where  f is  the  body  force,  r is  the  heat  source  and  we  are  using 
the  notation 


2%  = xt  + V-(xv).  (2.7) 

As  explained  in  the  introduction,  the  fields  (v(x,t), 

"t  (x, t)  , h(x, t) ) will  in  general  be  functions  of  bounded  variation, 
having  shock  discontinuities,  and,  consequently,  formulas  (2.3)-(2.6) 
will  only  be  satisfied  in  a generalized  sense  (in  the  sense  of 
measures  or  distributions) . In  particular,  the  classical  product 
differentiation  formula  does  not  generally  hold  for  functions  in  this 
class  so  that  (2.3)-(2.S)  cannot  be  simplified  in  the  standard 
fashion.  However,  in  the  special  case  of  a Lipschitz  continuous 
thermodynamic  process  (v(x, t) ,"t (x, t) , n(x, t) ) one  has  the  reduced 
form  of  the  balance  laws: 


+ yVP  = 

-"pV  .y 

(2.8) 

+ Vv'V  + 

Tvp  = f 

(2.9) 

+ v.Ve  + 

T p v-v  = r. 

(2.10) 

Since  x * 1/p,  one  obtains  from  (2.8) 


xt  + v-VT 


TV-v. 


(2.11) 


Furthermore,  using  (2.2)  one  easily  derives  from  (2.10) 
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n.  + v.Vn  = - (2.12) 

t ~ Q ' 

which  shows  that  the  Clausius-Duhem  inequality  (2.6)  holds  as  an 
equality  within  this  class  of  processes.  We  should  remark  here 
that  Lipschitz  continuous  processes  may  have  weak  (i.e.,  acceleration) 
waves  but  not  shock  waves. 


We  assume  that  (v(x , t)  , t (x , t)  ,"n(x,  t) ) is  a Lipschitz  con- 
tinuous process,  with  body  force  f(x,t)  and  heat  source  7(x,t), 
and  let  (v (x, t) , x (x , t) , n(x , t) ) be  any  other  admissible  process 
of  bounded  variation,  with  body  force  f(x,t)  and  heat  source  r(x,t). 
We  introduce  the  functions: 

H = P [j  (v-v) • (v-v)  + e - e + p(i-7)  - 0 (n-n)J  , (3.1) 

F = Hv  + (p-p)(v-v).  (3.2) 

We  will  be  using  H as  an  estimate  of  the  "distance"  between  the 
two  processes  while  F is  the  "flux"  of  H.  In  order  to  determine 
how  H evolves  in  time,  we  have  to  estimate  the  expression 
Ht  + V*F.  On  account  of  (2.7), 

Ht  + V*~  = I?  + v* t (P-P) (v-v)l . (3.3) 

Using  the  identity 

= * $ + X(¥t  + v-V*),  (3.4) 

which  holds  in  the  sense  of  measures  whenever  x is  any  function 
of  bounded  variation  and  41  is  Lipschitz  continuous,  one  obtains 
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6 i _ <5(PV)  _ 

Ht  + v’?  = + pe]  ' Y*”St PylYt  + VY'Y)  (3.5) 

+ py- tyt  + vy-v]  - p [£”t  + y-Ve]  + p + v • (pV) 

- Pp[Tt  + V-VT]  - pT[pt  + v-Vp]  - 0 ^-Pn^ 

- Pn[6  + v*V0]  + prf[0  + v-V0]  + p6  [rf  + yVrj] 

t t ^ t 

+ V-(pv)  - V-(pv)  - V • (pv)  + V-(pv). 

With  the  help  of  the  balance  laws  (2.3) - (2.6)  and  (2.9)  we  get 

Ht  + - -V- (pv)  + Pf -v  + pr  - v- [-Vp  + pf]  (3.6) 

- pv-[-7vp  + 7]  + pv-  [-TVp  + 7] 

- p(v  - y)-Vy (v  - v)  + pfc  + V • ( pv ) - pT[pt  + v-Vp] 

- pfe^  + v • v"e  + p[Tfc  + v • v"t]  - "nte  + v • V0 ] } 

■ Pr  g-  - P(n  - n)  [0  + v-vS"] 

♦ V-(pv)  - V-(pv)  - V-(pv)  + V-(pv). 

We  simplify  the  right-hand  side  of  (3.6),  noting,  in  particular, 
that,  by  virtue  of  (2.2), 
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— - ~ - ~ ~ ^-i.-sa:^^=3a:a 


thus  arriving  at 

H-t  + V*F  < p(f  - T)  • (v  - v)  - P(v  - y)*Vy-(v  - v)  (3.9) 

+ Pr  - Pr  g-  + (1  - pT)  [pt  + v-Vp]  + (p  - p)V-v 

- p(n  - "n)  [9"t  + v*V0]  - p(n  - n)  (v  - v)-ve 

= P(f  - £) • (v  - v)  - P(v  - Y)’7Y'(Y  ' Y) 

+ Pr  - Pr  ^ P(t  - T){|£—  [Tt  + v-VT]  + |£-  [nt  + v-Vn]} 

*\  a ^ ^ 

+ (p  - p) V • y - p(n  - [Tt  + y.vT]  + [nt  + v-Vn]} 

- p (n  - n)  (v  - v) • V0 . 

On  account  of  (2.2), 

3p*  _ 30*  _ 32e*  (3.10) 

3n  3t  " aian 

so  that,  with  the  help  of  (2.11)  and  (2.12),  we  may  rewrite  (3.9) 
in  the  form 
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Ht  + V*F  < P(f  - £ ) • Cv  - v)  - P(v  - v) -Vv- Cv  - v)  (3.11) 

■ tp  ■ p - (T  ■ T)  - (n  - n) ] v • v 

- p If  (T  - 7)Zv-v-  p (T  - Y)Cn  - n)v*v 

♦ p|  19  - f ^ (t  - T)  - |il  (i,  - n)! 

0 

. p X I (6  . e)2  + p I_i_E  (0  - ?)  - p(n  - n)(v  - v)-V0. 

F y 

Estimate  (3.11)  will  be  the  tool  for  establishing  stability. 
The  crucial  observation  is  that,  by  virtue  of  (2.2),  H and  F are 
of  quadratic  order  in  the  differences  (v  - v , t - x , n - n)  and 
that  the  right-hand  side  of  (3.11)  is  of  quadratic  order  in 
(v  - v,  T - T,  n - n,  f - T,  r - 7).  In  addition  to  that  we  have 
to  make  certain  that  H is  positive  definite.  To  this  end,  we 
impose  on  the  constitutive  equations  the  following  restrictions: 

cH>„  * »■  <!£>p  > »■ 

In  (3.12)  we  are  using  the  notation  of  classical  thermodynamics, 
visualizing  p as  a function  of  (i,n)  and  taking  the  partial 
derivative  with  respect  to  t and  then  visualizing  9 as  a 

function  of  (p,n)  and  taking  the  partial  derivative  with  respect 
to  n.  Assumption  (3.12)  is  physically  reasonable,  for  as  long  as 
the  fluid  does  not  undergo  any  phase  transitions,  and  its  connection 
with  Gibb’s  stability  has  been  established  in  classical  thermo- 
dynamics. In  particular,  (3.12)  is  satisfied  by  the  constitutive 
equations  of  polytropic  gases. 
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In  order  to  see  the  implications  of  (3.12)  on  the  function 
e*(T,n),  note  first  that,  by  virtue  of  (2.2), 


92£* 


(3.13) 


Furthermore,  by  chain  rule, 


l3n'\ 


rii) 

l3nJx 


3nJ 


8P  n 

■ 

■{- 


32e*  -1  a2e*  32  * 
3T  3i 


3n 


2*o 

(i_ "} 

^ 3 T 3 ly  ’• 


(3.14) 


It  is  clear  from  (3.13)  and  (3.14)  that  (3.12)  holds  if  and  only 
if  the  function  e*(x,n)  is  (locally)  uniformly  convex. 

Let  us  now  assume  that  the  fluid  is  confined  in  a vessel  with 
rigid  boundary  that  occupies  a smooth  bounded  domain  r c Rn  (n  = 1 
2 or  3).  This  leads  to  boundary  conditions 


Y*v  = 0 (3.15) 

where  v is  the  unit  normal  on  the  boundary  of  R . Our  stability 
result  is  given  in  the  following 

Theorem.  Let  (v(x,t),  T(x,t),  n(x,t))  be  a bounded  and  uniformly 
Lipschitz  continuous  process  defined  for  x e r,  t > 0 and 
satisfying  the  balance  laws  (2.3)-(2.10)  and  the  boundary  condition 
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v*  v = o.  Consider  any  admissible  process  (v(x,t),  x(x,t),  n(x,t)) 
of  bounded  variation  defined  for  x e fi,  t > 0 and  satisfying  the 
balance  laws  (2.3) -(2.6),  the  boundary  condition  v-v  = 0 and  the 
bounds 

| T (x,t) | , |p(x,t)|,  | n(x,t) | , |0_1(x,t) | < C,  x e fl,  t > 0.  (3.16) 


Then  there  are  positive  constants  A and  a,  depending  solely  upon 
bounds  of  (v,x,n)  and  its  Lipschitz  constant,  the  function  e*(x,n) 
and  the  constant  C in  (3.16),  such  that,  for  any  t > 0, 


(|v(x,t)  - v(x,t)|2  + |x(x,t)  - x(x,t)|2  + |n(x,t)  - n(x,t) |2}dx  (3.17) 

< Aeat  f ( | v(x, 0)  - v(x,0) | 2 + |t(x,0)  - ^(x,0)|2  + |n(x,0)  - n(x,0)|2}dx 

Ip  ~ 

+ ( (|f(x,s)  - r(x,s)|2  + |r(x,s)  - r(x,s) | 2}dxds  . 

J0  Jn 


Proof . Since  the  right-hand  side  of  (3.11)  is  of  quadratic  order 
(Y  - Y,  t - t,  n - n,  f - f,  r - 7),  we  have 

Ht  + v'5  < t»{  | v - v|2  + |x  - x|2  + |n  - n| 2 + |f  - f|2  + | r - r|2}  (3.18) 

with  b depending  solely  upon  bounds  of  (v,7,n)  and  its  Lipschitz 
constant,  the  function  e*(x>n)  and  C in  (3.16).  Integrating 
(3.18)  over  n x [0,t],  applying  Green's  theorem  and  noting  that  by 
(3.2)  and  the  boundary  conditions  we  have  F>v  = 0,  we  obtain 


Now  H is  of  quadratic  order  in  (v  - v,  t - T , n - o) 

so  that 


I H(x,0)dx  < M { | v (x , 0)  - v (x , 0) | 2 (3.20) 

Jn  Jn 

+ |t(x,0)  - T(x,0) | 2 + | n(x, 0)  - n(x,0) | 2}dx. 

On  the  other  hand,  since  £ (t,n)  is  uniformly  convex,  H is 
positive  definite  and  so 

' 2 

H(x,t)dx  > m ( | v (x , t ) - v (x , t ) j (3.21) 

Jn  Jn  " ~ 

♦ |T(x,t)  - T(x , t) | 2 + | n(x, t)  - n(x,t)|2}dx. 

The  positive  constants  M and  m in  (3.20)  and  (3.21)  depend  on  the 
same  factors  as  b in  (3.17).  We  now  combine  (3.19),  (3.20),  (3.21) 
and  we  apply  Gronwall’s  inequality  (e.g.  [10,  1.6.6])  thus  arriving 
at  (3.17)  with  u = b/m  and  A = max{M/m,b/m) . This  completes  the 
proof. 

The  preceeding  proposition  yields  immediately  the  following 
result  on  uniqueness: 
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Corollary.  Let  (v(x,t),  x’fx.t),  "rT ( x , t ) ) be  a uniformly  Lipschitz 
continuous  process  defined  for  x e ft,  t > 0 and  satisfying  the 
balance  laws  (2.8) -(2.10)  and  the  boundary  condition  v-v  = o.  Then 
there  is  no  other  admissible  process  (even  in  the  broader  class  of 
functions  of  bounded  variation)  that  satisfies  the  same  initial  and 
boundary  conditions  and  has  the  same  body  force  and  heat  source  as 
the  process  (v(x,t),  T(x,t),  n(x,t)). 

One  may  also  establish  uniqueness  and  stability  results  for  the 
case  where  the  fluid  occupies  the  entire  space  by  applying  the 
techniques  developed  in  [6] . 
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